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ABSTRACT 

We examine the effect of a threshold bias on the power spectrum and the bispec- 
trum in an ensemble of numerical simulations (Gaussian initial perturbations with 
power law spectra P(k) ~ k n } n = +1, 0, —1, —2) and compare our results with the- 
oretical predictions. Our simulations are evolved sufficiently that on the scale where 
we apply the threshold the rms fluctuation has developed significantly into the non- 
linear regime. Thus, predictions based on perturbation theory do not necessarily 
apply. Nevertheless, we find our results for the power spectrum, biased power simply 
amplified by a numerical factor, follow predicted trends, far beyond the regime where 
perturbation theory is expected to be valid. 

We find that the biased bispectrum continues to follow the so-called hierarchical 
form, with reduced three-point amplitude Q 1 in the strongly nonlinear regime, 
independent of initial spectrum. In the quasi-linear perturbative regime the three- 
point amplitude depends on configuration shape, a behavior that is found to give 
useful information about the amount of bias without information about the unbiased 
matter distribution. 

Subject headings: methods: numerical 
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1. Introduction 

Soon after Kaiser (1984) showed that the statistical properties of Abell clusters could be 
at least qualitatively understood as properties of high density peaks in the distribution of 
galaxies, it was recognized that the galaxy distribution itself might be similarly biased relative 
to the underlying distribution of mass. A bias in galaxy formation makes the statistical 
properties of the galaxy distribution different from those of the underlying mass distribution. 
In the intervening years a large number of statistical and physical bias mechanisms have been 
proposed, which usually but not always lead to a galaxy distribution more strongly clustered 
than the matter distribution. Results can be summarized by the bias parameter 6, where 
the correlation function of the biased galaxy distribution ^ is related to that of the matter 

distribution £ by ^ = b £. It is almost axiomatic that some bias is necessary to reconcile the 
smoothness of the microwave background and the observed clustering strength of galaxies. 

Theoretical calculations of what to expect from bias in galaxy formation are confined 
mainly to the perturbative regime. To study the further implications when fluctuation am- 
plitudes evolve into the nonlinear regime the standard technique is numerical simulation. 
Melott & Shandarin (1993) have investigated the power spectrum and Fry, Melott, & Shan- 
darin (1993) the bispectrum of the mass distribution at various evolution stages for an 
ensemble of three-dimensional numerical simulations for a series of scale-free initial condi- 
tions. In this paper we study the power spectrum and bispectrum after applying a threshold 
bias to the same ensemble of models (Gaussian initial perturbations with power law spectra 
P(k) ~ k n } n = +1, 0, —1, —2). In Section 2 we review the behavior of the power spectrum 
and the bispectrum expected in perturbation theory, including effects of gravity and non- 
linear bias. Section 3 presents our results from numerical simulations. Section 4 contains a 
discussion of what we have learned. 



2. Theoretical Predictions 
2.1 Gravitational Instability 

The large-scale distribution of matter in the universe is shaped by the action of gravity 
on small initial seed fluctuations, usually taken to be Gaussian. Gravitational instability is 
nonlinear, thus precluding a completely analytic treatment. However, on sufficiently large 
scales, density fluctuations are weak enough that perturbation theory often suffices. To 
linear order, the fluctuation amplitude <5(fc,t), the Fourier transform of the density contrast 

8 = [p(x } t) — p]/p } grows by an overall scale factor, 8 = A(t)<5g(fc), where <*>o(fc) is the 
amplitude at some early time to (see Peebles 1980, § 10-13). For the canonical model, matter 

dominated, ft = l } A = k = p = = 0, this factor is Ait) ~ t^l^ ~ a(t), where ait) is the 
cosmological expansion factor, and normalization can be chosen so that Ait) = ait). Thus, 
perturbation theory predicts to leading order the power spectrum P(k) = ( |<5(fc)|^ ) simply 
grows with time, 

P(k) = a 2 (t)P (k). (1) 

For Gaussian initial conditions, all higher order irreducible correlations for 8 vanish ini- 
tially. However, in a nonlinear theory, even for a Gaussian initial distribution, nonlinearities 
induce nonvanishing higher order correlations for all orders (see Fry 1984, Goroff et al. 1986, 

Bernadeau 1992). For the three-point function in particular, -B123 = ( 8(ki)8(k2)8(k^) ) for 
^2 fcj- = 0, gravitational instability gives 



#123 = #12 + #13 + #23 



(2) 
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where 

10 t f (hi kj\ 4 



Pik^Pih). (3) 



Note that since P ~ a?(t), B ~ cfi. 

We can remove the main dependence on scale, on time, and on the initial power spectrum 
by using the hierarchical amplitude Q, defined as 

Qm = PlP2 + P1P3 + P2P3' W 

In a pure hierarchical model, Q would be strictly constant. For equilateral triangle config- 
urations kj/kj = 1 and fcj • kj = —1/2 for all pairs, and equation (2) gives Ql23 = 4/7, 
independent of k, independent of time, and independent of the initial power spectrum. For 
a power law initial P(k), Q123 is independent of time and of overall scale k, but depends 

on configuration shape in the ratios k^/kj and angles fcj • kj. This dependence of the three- 
point amplitude on configuration shape is the signature of structure formed by gravitational 
instability. 

2.2. Bias 

Over the past decade it has become popular to allow for a bias in galaxy formation. 
Kaiser (1984) introduced the concept of a threshold bias originally to relate the clustering 
properties of Abell clusters to those of galaxies by examining the statistical properties of those 
regions where the density contrast exceeds some threshold, 8 > t. Expressing the threshold 
as a multiple of the rms fluctuation, t = vcfq (possibly after smoothing the distribution), 
Kaiser showed that, for v > 1 and £q < 1, the biased distribution has a correlation function 
Cb = ^;/C/Co? with b v = v. This picture was soon extended to include the possibility that the 
galaxy distribution is biased relative to the mass density. The implications of a threshold bias 
have been investigated analytically in great detail for an underlying Gaussian distribution 
(Kaiser 1984; Politzer & Wise 1984; Bardeen et al. 1986; Jensen & Szalay 1986). Jensen & 
Szalay extended the results to small z/, obtaining the amplification factor 

K = -^- y F(u) = y§ * e" 2 / 2 edc(u/V2). (5) 

The denominator F(v) approaches 1 when v is large, thus reducing to the Kaiser result, 
and vanishes linearly as v — > 0, such that the factor multiplying <£/£o is = 2/7T. Many 
other more or less physically well-motivated models have been considered, with results often 
summarized by the value of a bias parameter 6, the ratio of the density contrast in the biased 
distribution to that in the underlying density, 8^ = bS. If a linear bias were the complete 
story, the spectra would obey P^ = b^P, = b^B, and thus = Q/b. However, it seems 
likely that the process of galaxy formation leads to a galaxy distribution which does not 
depend linearly on the underlying mass distribution 

The threshold model is one example of a nonlinear bias. Applied to an underlying 
Gaussian, the threshold model gives Q = 1 in the weak fluctuation regime (Politzer & 
Wise 1984; Jensen & Szalay 1986), while for strong clustering Q — » 1/3 (Jensen & Szalay 
1986). Generally we would expect that the galaxy distribution is a functional of the matter 
density, 8j ) (x) = f[8(x)]. If the range of effects that determine the efficiency of galaxy 
formation is not too large, this perhaps can be approximated as an arbitrary but local 
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function, 6^ = f(S). Szalay (1988) has considered an arbitrary biasing function, expanded in 
Hermite polynomials as /(<*>) = ^2 Jj e Hj e (6) / k\ } applied to a Gaussian underlying density field 
and obtained a hierarchical three-point function with Q = J^/ jf. The next step is to account 
for a non-Gaussian matter distribution, as that produced by gravitational instability. Fry & 
Gaztanaga (1993) considered the implications of an arbitrary local bias function applied to 
an underlying hierarchical mass distribution. Expanding /(<*>) in a series 

= £ yhs k , (6) 

k 

they found, to leading order for small fluctuations, hierarchical galaxy correlations for all 
orders, with amplitudes that depend on the original Q n and on the bias parameters b n . In 
particular, the biased three-point function is hierarchical, with amplitude 

1 &2 

Qm,b = jQm + ( 7 ) 

(Fry & Gaztanaga 1993; Juszkiewicz et al. 1993), showing a contribution from nonlinear 
gravity that depends on configuration and a contribution from nonlinear biasing that is 
independent of configuration. No model of bias that has been suggested to date generates a 
contribution to Q that depends on configuration, although a more complicated bias model 
could conceivably do so. Further dependence on configuration shape can arise from higher 
order moments for non-Gaussian initial conditions (Fry & Scherrer 1994). 



3. Numerical Results 
3.1. The Power Spectrum P(k) 

Melott & Shandarin (1993) and Fry, Melott, & Shandarin (1993) have investigated the 
power spectrum and the bispectrum of the mass distribution for an ensemble of numer- 
ical simulations. The results in the perturbative regime agree within uncertainties with 
gravitational instability predictions in perturbation theory and obey expected scalings in 
the nonlinear regime up to departures that we believe are understood. In this paper we 
examine the effects of a threshold bias applied to the same ensemble of models. The simula- 
tions are performed on a staggered 128^ mesh (Melott 1986), starting from small amplitude 
Gaussian initial conditions with random phases for the Fourier components. For the differ- 
encing scheme used, growth rates of fluctuations are suppressed for wavenumbers larger than 
k > k-^y/2 = 32 (considerably better than the usual non-staggered mesh value ~ &N y /4). 
We present results for wavenumbers in the range 1 < k < 32, where k = 1 corresponds to 
the fundamental mode. 

Models are identified by initial power spectrum index, Pq ~ k n } for n = +1, 0, — 1, 
and —2 out to the Nyquist frequency k-^y = 64. We have also performed simulations for a 
spectrum with n = —3, but as our earlier results bear out, results of any attempt to realize 
such a divergent spectrum on a finite volume are severely affected by the cutoff, and so in 
this paper we examine only n > —2. For each initial spectrum we have information for a 
succession of evolution stages labeled by k n \, the scale of nonlinearity given by 

nl d 3 kP(k) = 1, (8) 
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evaluated for the analytically evolved initial spectrum, P = a (t)PQ(k). A given k^\ repre- 
sents evolution for an expansion factor of approximately a = Aq~^ (&Ny / k n i)^~^ n ^ 2 , where 
Aq is the initial rms fluctuation in one cell. 

For each initial spectrum and for each stage of evolution we construct a "biased" dis- 
tribution by keeping only particles in cells with density above an absolute threshold 8 > t 
on the scale of f cell at three different threshold levels, t = 3, t = fO, and t = 30. Table f 
provides the rms cell fluctuation erg at each stage necessary to translate this threshold into 
a number of standard deviations t = v(tq. At evolution stage k n \ = 32 fluctuations on the 
scale of f cell are just becoming nonlinear (ctq = 1.2 -1.8), while at stage k n \ = 4 the degree 
of nonlinearity is large (ctq = 8.4-12.3). 

After applying the threshold, unlike the situation in our earlier analyses of the unbi- 
ased distribution, the number of particles remaining is so small that we now must remove 
discreteness contributions to the the Fourier spectra, 

(\6(k)\ 2 )=P(k) + ± (9) 

( hhh ) = £123 + jj[P(h) + P(h) + P(h)] + ^ , (io) 

where P and B are the contributions from clustering. The numerical data suggest the 
appropriate discreteness measure is the number N = N c of cells occupied above threshold. 
For comparison with later results, Table 2 shows the discreteness correction Pj = 1/N C 
averaged over two random realizations for k n \ = 32 and k n \ = 4 (the earliest and latest 
evolution stages), as a function of spectra index n and threshold t. For late stages of evolution 
at lower thresholds, N is large, and the discreteness correction is insignificant. At early stages 
of evolution and higher thresholds, the correction can be large. For the most extreme case, 
n = +1, k n \ = 32, t = 30, we have for two random realizations N c = 3 and N c = 0. It is 
likely in this case that the excess power above discreteness is just random noise; we include 
it because there may be information in the distribution of this noise. 

Figures 1-4 show P(k) after the correction for discreteness. Error flags in the figures 
show the dispersion between random realizations, except for k n \ = 32, t = 30, where the error 
indicates the variation within the one realization that has N c ^ 0. Recall that early stages 
with large thresholds are dominated by discreteness and are not likely to be reliable. Figures 
1 and 2 both show results for initial spectrum n = — 1. Within each window in Figure 1 the 
four curves show the effect of increasing threshold at the same evolution stage, as indicated. 
At an early evolution stage (k n \ = 32) the power depends strongly on the threshold; however 
his dependence steadily diminishes with evolution, becoming almost negligible at k n \ = 4. 
Qualitatively, this behavior could be anticipated: at later stages a considerable amount of 
mass is in dense clumps, and increasing the threshold makes almost no difference in the 
mass distribution. The four windows in Figure 2 show increasing threshold, from no bias to 
t = 30, while the curves within each window show stages of evolution. The effect mentioned 
above is demonstrated in a different form in Figure 2, where the effect of time evolution 
steadily diminishes with increasing threshold. That the short wave part of the spectrum 
of the highest peaks (8 > 30) appears to decrease with time is probably an artifact of the 
shot noise subtraction, which dominates the power here. Figures 3 and 4 show the effect of 
increasing threshold for the four values of initial spectral index, for k n \ = 32 (Fig. 3) and 
at k n \ = 4 (Fig. 4). The dramatic difference between Figure 3 and Figure 4 illustrates the 
influence of the threshold bias in the mildly nonlinear regime (erg = 1.2-1.8 in Fig. 3) and 
highly nonlinear regime (ctq = 8.4-12.3 in Fig. 4). 

These results follow the qualitative trend that a larger bias threshold produces a greater 
amplification of the power spectrum. Surprisingly, the linear bias prediction P\){k) = b^P(k), 
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is at least approximately quantitatively valid as well. The first column under each heading 
in Table 3 shows the value of b derived from the power spectrum averaged for 1 < k < 32 
(this is not necessarily a good fit) and the second the Jensen & Szalay (1986) predicted 
amplification factor b v for the given threshold (eq. [5]). Comparing the two columns one 
can see a general correspondence, although a difference of a factor of 2 or greater is not 
unusual. Note that this agreement obtains only if we omit the factor erg that appears in the 
denominator of the perturbation theory prediction. 

3.2. The Three-Point Amplitude Ql23 

We next present the reduced three-point amplitude <5l23- Figures 5-8 show Q(k) for 
equilateral triangle configurations after the discreteness correction in a manner similar to 
the power spectrum. For clarity, errors are shown only for selected curves; in general, 
errors are comparable to the apparent fluctuation from point to point. As described in § 2, 
the gravitational instability prediction for the quasi-linear regime is Q = 4/7, indicated as 
the horizontal dashed lines on the left of each panel; and numerical simulations suggest 
Q(n) 3/(3 + n) for no bias in the nonlinear regime, indicated by the horizontal dashed 
lines on the right. As in Figure 1, the four windows in Figure 5 show increasing evolution for 
n = — 1, from k n \ = 32 to k n \ = 4. The curves in each window show increasing threshold, as 
labeled in the caption. In the quasi-linear regime k < k n \ } Q 4/7 appears to hold within 
the errors (which are large for small k). ^From equation (7), this is expected for no bias but 

also holds for certain combinations of bias parameters; for instance, for b ^> 1 if b^/b^ 4/7. 
For k > k n \ } Q steadily grows, towards the expected Q = 3/2 with no bias, but with a mild 
but consistent decrease with increasing bias, closer to Q = 1 for t = 30. The four windows 
in Figure 6 show n = — 1 with increasing threshold, from no bias to 6 > 30; the curves in 
each window show increasing evolution, as labeled in the caption. With no bias the rise from 
Q = 4/7 towards Q = 3/2 is noticeable, especially at late evolution, while with the largest 
bias t = 30, Q = 1 for scales in the nonlinear regime. 

The windows in Figures 7 and 8 show increasing bias for different values of initial spectral 
index, from n = +1 to n = —2, at k n \ = 32 (Fig. 7) and at k n \ = 4 (Fig. 8). Fry, Melott, 

& Shandarin (1993) show that errors in Q go as (k n \/ k)^~^~ n \ and thus are very large for 
positive n and small k; points with AQ ^> Q are not shown. Figure 7 shows an early stage 
when all scales shown are in the quasi-linear or mildly nonlinear regime. Therefore we do 
not see much departure from Q = 4/7, except for n = —2 with no bias. In Figure 8, on the 
other hand, all scales except the smallest k are in the strongly nonlinear regime, and we see 
a steady growth of Q from Q = 4/7 for small k towards Q(n) = 3/(3 + n) at large k without 
bias or Q 1 at high threshold. The models dominated by the small scale perturbations, 
n = +1 and n = 0, do not display much dependence on the bias factor; however, in the 
n = — 2 model this dependence is quite significant. Figure 9 shows one final view of different 
parameter combinations, four different values of the initial spectral index n in windows of 
increasing bias threshold. As shown in Figure 9, with no bias Q approaches the corresponding 
value Q(n) = n/(n + 3) at k = 32; with the largest bias, 6 > 30, all curves look similar and 
approach the value Q = 1 at k = 32. The no bias window also shows k n \ = 8. 

Finally, we examine the dependence of Ql23 on configuration shape. We plot in Figure 10 
the reduced three-point amplitude Q(6) for configurations with k\ = 32 and k^ = 16, 
separated by angle 6 } for initial n = — 1. For clarity error flags are shown only for no bias 
and for t = 3. The earliest stage shows a substantial variation of Q with 6 } with enhanced 
correlation for 6 = and 6 = tt. One can intuitively understand the origin of such a shape 
dependence from gravitational instability. The Zel'dovich (1970) approximation describes 
pancaking (highly anisotropic collapse) from initial conditions with negligible power at large 
k in the initial conditions. Coles, Melott, & Shandarin (1993; see also Melott, Pellman, & 
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Shandarin 1994) showed that a simple extension ol the Zel'dovich approximation applies to 
Gaussian initial conditions with a wide variety of power spectra. Thus the enhancement of 
local anisotropy is a generic result of gravitational instability. In transform space, the effect of 
a two-dimensional spatial density enhancement is a preference for wave vectors orthogonal to 
the two-dimensional structure, and in statistics averaged over orientations this survives as an 
enhancement for nearly colinear configurations. The dumpiness of pancakes and filaments 
eliminates this effect at late stages. Since the size of the triangle in the transform space 
(k\ = 32 and k^ = 16) remains constant during the evolution, we probe scales a little greater 
than the size of clumps at the stages k n \ = 32 and 16 and see enhancement of Q at small 6 
for n = —2. On the other hand at later stages k n \ = 4 the triangle probes the interior of the 
clumps and the enhancement disappears. 



4. Discussion 

In this paper we have examined the power spectrum and bispectrum for an ensemble of 
numerical simulations with an applied threshold bias. As predicted in perturbation theory we 
have seen that that the general effect of a threshold bias on the power spectrum is to amplify 
the power by a roughly constant factor that can be written as b ; this result applies from 
the quasi-linear regime into the strongly nonlinear regime. The numerical results follow 
qualitatively the trend predicted in the amplification factor b v of equation (5) (Jensen & 
Szalay 1986), but depart from exact numerical agreement by factors of 2 or larger. A more 
serious departure from the prediction is that agreement is achieved only when the factor £q 
in the prediction £g = &^£/£o i s se ^ to 1. Yet it is perhaps the agreement rather than the 
departure that is the more surprising, because many of the numerical results are far outside 
the domain where perturbation theory should apply. 

The effect of bias on the three-point function is different in the quasi-linear and strongly 
nonlinear regimes. In the nonlinear regime the effect of bias on the amplitude Q(k) for 
equilateral triangles is to modify the value Q(n) = 3/(3+n) for no bias to Q 1, independent 
of n, for high bias at large k. This confirms the result that biased Q 1 found previously 
by Melott & Fry (1986) for simulations on a much coarser grid. This value is expected for 

the threshold bias when b ^> 1, since effectively b^/b 2 = 1, but it appears to hold for large 
threshold even when the amplification b is not much different from 1. In the quasi-linear 
regime, where we ought to be able to make reliable predictions, the errors become large and 
it is difficult to draw specific conclusions. It seems that Q(k) again approaches a constant 
value for small k, but it would be difficult to conclude that this value is different from the 
gravitational instability result Q = 4/7. Note that even so, for this to be consistent with 
equation (7) with a large b inferred from the power spectrum P(k), a b^ term is necessary. 

Further information can be extracted from the dependence of the three-point amplitude 
Q on configuration shape. Results for k n \ = 32 are the least evolved and thus should be the 
closest to perturbation theory predictions. The perturbation theory prediction depends on 
spectral index, which evolves slightly from the initial value n = — 1 and depends slightly on 
scale: for small wavenumbers, 1 < k < 16, n = —1.2; while for 16 < k < 32, n = —1.6. 
For t = 3, fitting Q(6) to the form in equation (7) with n = —1 gives b = 4.5 ± 0.2, 
b 2 /b 2 = 0.51 ± 0.01. For n = -1.5, a value we might infer from P(k) at the stage observed, 
b = 5.5 ± 0.3, b 2 /b 2 = 0.52 ± 0.01. For t > 10, a fit with n = -1 gives b = 6.9 ± 0.6, 

b 2 /b 2 = 0.67 ± 0.01; n = -1.5 gives b = 8.1 ± 0.8 and b 2 /b 2 = 0.68 ± 0.01. (These errors are 
formal errors in a linear least squares fit, ignoring parameter correlations.) For comparison, 
the ratio of the biased and unbiased power spectrum in these two cases over the same range 
in k gives b = 4.10 ± 0.01 and b = 9.5 ± 0.2. It appears that the shape dependence of 
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the three-point function slightly overestimates the amount of bias, but we note that even 
for k n \ 32, the results have evolved beyond the truly perturbative regime. In the highly 
evolved results, k n \ = 4, there is less dependence of Q on configuration shape, even with no 
bias. What dependence Q has on 6 is very different from that in perturbation theory, but 
it appears to be robust and to be affected by bias only in amplitude, with no offset [&2 = 
in eq. (7)]. Using the form = Q/b applied to the curves for k n \ = 4 in Figure 10 gives 
b = 1.24 for t = 3, b = 1.36 for t = 10, and b = 1.54 for t = 30, values that compare well 
with those obtained from the power spectrum in Table 3. 

In conclusion, we have seen that a threshold bias, even at a very high level, changes the 
power spectrum and bispectrum of the biased distribution, but in ways that make it still 
possible to extract useful information about the underlying matter distribution and about 
initial conditions. Bias does not change the shape of the power spectrum except when the 
threshold is very large. In the nonlinear regime, with almost any appreciable bias Q I 
for equilateral triangles, independent of initial spectrum or degree of evolution. However, 
in the perturbative regime, k < k n \ } the characteristic dependence of the bispectrum on the 
shape of the triangular configuration in perturbation theory (eqs. [2-3]) allows in principle 
the separation of the contributions of bias and gravity, and provides an important new, 
totally independent method for estimating b that does not require separate knowledge of the 
amplitude of the unbiased power spectrum. 
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Cambridge University, and JNF to Fermilab for their hospitality during the period in which 
some of this work was done. We thank the National Center for Supercomputing Applications, 
Urbana, Illinois for the grant of Convex C3 time which made this work possible. 
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Table 1 
R.M.S. Cell Fluctuation erg 



n k n \ = 32 k n \ = 16 k n \ = 8 k n \ = 4 



+ 1 1.20 3.06 6.42 12.30 

1.22 2.95 5.91 10.74 

-1 1.22 2.67 5.83 9.83 

-2 1.84 3.19 5.16 8.36 



Table 2 

Discreteness Power Pj = 1/N ( 



k n \ t n = +1 n = n = — 1 n = —2 



32 3 1.65 x 10" 5 1.33 x 10" 5 1.92 x 10" 5 1.50 x 10" 5 

10 5.02 x 10" 4 2.20 x 10" 4 2.95 x 10" 4 1.14 x 10" 4 

30 3.33 x 10" 1 4.63 x 10" 2 1.25 x 10" 2 1.31 x 10" 3 

4 3 1.86 x 10" 5 1.60 x 10" 5 1.52 x 10" 5 1.26 x 10" 5 

10 3.26 x 10" 5 3.03 x 10" 5 3.32 x 10" 5 3.80 x 10" 5 

30 6.98 x 10" 5 7.15 x 10" 5 8.48 x 10" 5 1.19 x 10" 4 



Table 3 

Bias Amplification Factor b (n = — 1) 



t k n \ = 32 k n \ = 16 k n \ = 8 k n \ = 4 



3 4.1" 2.9 & 2.1 1.6 1.4 1.2 1.2 1.0 

10 9.1 8.3 3.2 4.0 1.8 2.1 1.5 1.5 

30 24.1 24.6 5.2 11.3 2.5 5.8 1.7 3.3 



'Observed b from P(k) 
Predicted b v from equation (5) 
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. 1. — P(k) vs. k for models with n = —1. The four windows show evolution stages 
k n \ = 32, k n \ = 16, k n \ = 8, and k n \ = 4. In each window curves are plotted for no bias 
(solid line), and for bias thresholds 8 > 3 (short-dash line), 8 > 10 (long-dash line), and 
8 > 30 (dot-dash line). 
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Fig. 2. — P(k) vs. k for models with n = —1. The four windows show no bias and bias 
thresholds 6 > 3, 6 > 10, and 6 > 30. In each window curves are plotted for evolution 
stages k n \ = 32 (solid line), k n \ = 16 (short-dash line), k n \ = 8 (long-dash line), and 
k n \ = 4 (dot-dash line). 
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Fig. 3. — P(k) vs. k for models at evolution stage k n \ = 32. The four windows show initial 
power spectrum index n = +1, n = 0, n = — 1, and n = —2. In each window curves are 
plotted for no bias and for bias thresholds 6 > 3, 6 > 10, and 6 > 30, indicated as in 
Figure 1. 



14 




. 4. — P(k) vs. k for models at evolution stage k n \ = 4. The four windows show initial 
power spectrum index n = +1, n = 0, n = — 1, and n = —2. In each window curves are 
plotted for no bias and for bias thresholds 6 > 3, 6 > 10, and <5 > 30, as in Figure 1. 
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Fig. 5. — Q(k) vs. k for models with n = —1. The four windows show evolution stages 
k n \ = 32, k n \ = 16, k n \ = 8, and k n \ = 4. In each window curves are plotted for no 
bias (solid) and for bias thresholds 8 > 3 (short-dash), 8 > 10 (long-dash), and 8 > 30 
(dot-dash). Horizontal dashed lines show Q = 4/7 expected from perturbation theory 
and Q = 3/2 expected in the nonlinear regime with no bias. For clarity error bars are 
shown only on one curve in each panel. 
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. 6. — Q(k) vs. k for models with n = —1. The four windows show no bias and bias 
thresholds 8 > 3, 8 > 10, and 8 > 30. In each window curves are plotted for evolution 
stages k n \ = 32 (solid), k n \ = 16 (short-dash), k n \ = 8 (long-dash), and k n \ = 4 (dot- 
dash). 
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Fig. 7. — Q(k) vs. k for models at evolution stage k n \ = 32. The four windows show initial 
power spectrum index n = +1, n = 0, n = —1, and.n = —2. In each window curves are 
plotted for no bias and for bias thresholds 8 > 3, 8 > 10, and 8 > 30 as indicated in 
Figure 5. Horizontal dashed lines on the right-hand side of each panel show Q = 3/(n + 3) 
expected in the nonlinear regime with no bias. 
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. 8. — Q(k) vs. k for models at evolution stage k n \ = 4. The four windows show initial 
power spectrum index n = +1, n = 0, n = — 1, and n = —2. In each window curves are 
plotted for no bias and for bias thresholds 8 > 3, 8 > 10, and 8 > 30, as indicated in 
Figure 5. 
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Fig. 9. — Q(k) vs. k for models at evolution stage k n \ = 4. The four windows show no bias 
and bias thresholds 8 > 3, 8 > 10, and 8 > 30. In each window curves are plotted for 
n = +1 (solid line), n = (short-dash line), n = —1 (long-dash line), and n = —2 
(dot-dash line). The no bias window also shows k n \ = 8. Horizontal dashed lines on the 
right side of the no bias panel show the expected asymptotic value Q(n) = 3/(3 + n). In 
the remaining panels Q = 1 is marked. 
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Fig. 10. — Q(0) vs. 6 for models with n = —1, configurations with k\ = 32, k^ = 16, separated 

o o o 

by angle 9 such that = kf -\- + 2k\k^ cos The four windows show evolution stages 
k n \ = 32, k n \ = 16, k n \ = 8, and A; n j = 4. In each window curves are plotted for no bias 
(solid line) and for bias thresholds 8 > 3 (short-dash line), 8 > 10 (long-dash line), and 
<5 > 30 (dot-dash line). 



